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Abstract 

We introduce the chemical potential in a system of two-dimensional 
massless fermions, confined to a finite region, by imposing twisted bound- 
ary conditions in the Euclidean time direction. We explore in this simple 
model the application of functional techniques which could be used in 
more complicated situations. 
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1 Introduction 



The thermodynamics of hadronic matter has recently received great attention 
[0, mainly in connection with the possible occurrence of a deconfining phase 
transition at finite temperature. The difficulty of studying this transition in the 
framework of QCD comes from the fact that confinement has not been deduced 
in this theoretical context. 

Consequently, effective models as, for example, the bag models |2[ ^ have 
been introduced to mimic the confining properties of strong interactions. In 
these models, fields are confined to bounded regions and subject to adequate 
boundary conditions. Finite volume effects turn out to be relevant, and they 
have been studied in the thermodynamical limit, for instance, in and for a 
nonvanishing chemical potential in . 

A complete analysis of the free energy for an MIT bag model at T > 0, 
with /i = 0, has been performed in 1^. The corrections introduced by chiral 
boundary conditions have also been studied 0. Functional methods have been 
used in these papers to isolate the finite temperature dependent pieces from the 
(divergent) Casimir energy. 

The aim of the present paper is to study the possibility of introducing a 
nonvanishing chemical potential through "twisted" boundary conditions in the 
Euclidean time direction and treat it following the methods developed in ^, ^ . 
As a first approach to this problem we present here the evaluation of the Gibbs 
free energy for the simple model of massless 1 + 1-dimensional fermions confined 
to a segment, for T and /i ^ 0. 

Even though, for this toy model, the eigenvalue problem can be exactly 
solved, we will rather follow an alternative approach. We will relate differences 
of free energies of the system to the Green function satisfying adequate (spa- 
tial and temporal) boundary conditions. This approach is more likely useful 
in the realistic four dimensional case, where eigenvalues cannot be explicitely 
solved for. In section 4, we will reobtain the results in Section 2, making use 
of the functional method developed in |^ , which is based on the evaluation of 
p-determinants. 



2 Two-dimensional fermions with /i 7^ 0. 

We consider a system of two-dimensional fermions confined in the segment [0, L] 
and subject to given boundary conditions. As is well known, the chemical 
potential can be introduced as an imaginary constant temporal gauge field 1^. 
If the temperature is T and the chemical potential is /i, the Grand canonical 
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partition function can be expressed as 

S(r,L,M) = e-^G(T,L,^) 

= / 2?^ efo '^*/o d^HDi0,L)-i^^J°)^ (1) 



where 



with 



D{P,L) = ^j°dt + ^j'd., for 0<t,x<l, (2) 



and "6c" means that we must cvahiatc the functional determinant of the differ- 
ential operator defined on functions satisfying the boundary conditions: 

B'4}{t,x) = 0, for X = 0, 1, 

(4) 

V'(l,a;) = -V(0,a;). 

Here B denotes an elliptic boundary condition to be satisfied at the spatial 
edges. 

So, we have the following relation for the Gibbs free energy: 



L, m) = rr{|^ In {D{(3, L) - i/.7°),J 
= Tr{-i^''Ki^{t,x;t',x')], 



(5) 



where Kbc{t,x;t',x') is the Green function of the problem, satisfying 
(£)(/?, L) - ifi-f") Kbc{t, X- 1', x') = 6{x - x')6{t - t') 

BKbc{t, x; t', x') = 0, for x = 0, 1, (6) 

Kbc{l,x;t',x') = -Kbc{0,x;t',x'). 
This can also be expressed in terms of 

D{f3, L) - = ef'f^^DiP, L) e""'^*, 

Kbc{t,x;t',x') = e^'f*k{t,x;t',x') e""''*'. 



(7) 
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where k{t,x;t' ,x') is the Green function of the problem 

D{(3, L) kit, x; t', x') = 6{x - x')6{t - t') 

Bei^l^^k{t, x; t', x') = 0, for a; = 0, 1, (8) 
A:(l, x; t', x') + e-f'l^kiO, x; t', x') = 0. 
This function has the development 

oo 

k{t,x;t',x') = L J2 kn{x,x') e-'^-'^^'-*'\ (9) 

n=—oo 

with the frequencies given by 



(2n+l)7r ^ 
uJn = , for n G Z. 



(10) 



The coefEcients satisfy 

kn{x, x') = e''^"-^^A;„(0, x') - ij^e-^'^"^^''-'''^H{x, x'), (11) 

with 

Now, wc must impose the boLmdary conditions on the spatial edges. We will 
adopt a static "bag-like " condition 

(B^) {t, x) = (1 + ii)ii){t, x) = 0, for a; = 0, 1, (13) 

such that, at a; = 0, 

(1 - 7')fcn(0, X')=(^\ 1 ) fcn(0, X') = 0. (14) 

This implies that 

fc„(0,x')= ( J ) ®( f{x') g{x') ). (15) 
On the other hand, at a; = 1, 

(l+7i)fc„(l,a;') = 0, (16) 
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which means that 



(1-0 

This gives 

oo 

k{t,x;t',x') = iL J2 



j^n^Lco M ® ( f{x') g{x') ) - j^ie-75Sini(i-x/) 



= 0. (17) 



2 cosh(J7„L) 

_^lg-75f2ni(a:-a;')//" (a; _ x')^ g-in„/3(t-f') _ 

So, we have for the derivative of the free energy 

-l3^{p,L,fx)=Tr{-ij^e''^*k{t,x;t',x') e-^^^t'} 

g— 750„L(1— X— x/) 



(18) 



.Hl3{t-t') 



E 



2 cosh(r2„L) 



+7^ 



— 75nTi-L(l+a:— a:/) 

2 cosh(il^^L) 



(19) 

For < x,x' < 1, the first term in the series converges absolutely and uniformly, 
even at {t, x) — {t', x'). So, it can be summed up in any order; in particular, one 
can first take the trace as a matrix, obtaining a vanishing result. This is not 
the case for the second term, containing the 7^-matrix. In fact. 



E 



E 

n=0 



+ E 



2 cosh(r2„X) 

-75n„L(l+a;-x/) 

2 cosh(r2„L) 



_ g-75n„i,(x-x/)^(2. _ J,/-) ] g- 



_ g— 760„I/(a:— X/) 




1 



-in„/3(t-t') _ 



-ia„/3(t-t') 



—')b^nL{l+X—Xl) 

2 cosh(f2„L) 



1 




-Ee- 

n=0 



-1b^nL{x — Xl) 



-H{x - x') 

H{x' - x) 



n=— 00 



L{x-xi) I ^{^' ~ ^) 

-H{x - x') 



^-ia„D{t-t') (20) 

g-iO„/3(t-t') 

-if2„/3(t-t') 



where the first two terms in the right hand side are absolutely convergent even 
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at x' = X. For the last two terms in this equation we have, for x' x, 



-■y^QnL^x—xf) 



-H{x-x') 
H{x' - x) 



-in„0{t-t') _ 



/ (iix-K/!3)[L{x-x')+il3{t-t')] 

' -H{x-x') 



2% 



^ _ ^--0-[L(x-x')L+i/3it-t')] 



(itJ.-7r/l3)[L{x' -x)L+if3{t-t')] 

H(x' - x)- 

\_ ^-'^[L(x'-x)L+iP(t-f)] J 
(21) 



and 



5-= E e 

n= — oo 



' H{x'-x) 



-(i/,+ ;r/ n'L(.r'-.r)-, .;(/-/' 



H{x' -x) 

-H{x-x') 



-iO„/3(t-t') 



27r 



V 



1-e 




[L(x'-x)-i/3(t-t')] 



-if (a; - x') 



-{ilJ.+w/l3)[L(x-x')-il3(t-t')] 



2Wr 



^ _ ^—0-[L{=^-^')-il}{t-t')] ) 

(22) 

which are singular at {t',x') = {t,x). Calling z = L{x — x') + i(5{t — t'), z = 
L{x — x') — i(3{t — t'), we have 



S+ + S- = 



p I - + + 7r//3 {H{x' - x) - H{x - x')) + 0{z/P^) ; ^ 
^ 1 

; -- in + Tr/ 13 {H{x' - x) - H{x - x')) + 0{z/p'^) 
z 
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(23) 



Now, we can evaluate 



L tr ^^o^^f3{t-t')^i + 5_)} = 

= _iM^/^/3(t-t') / 1 _ i + 2i/x + 0(| I //32) 
1^ 2: z 

^+0{\x-x'\Lyp) 



2tt 



(24) 
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This correspond to the foUowing (finite) contribution to ~(3^{T,L, ^) 



L 



= {t,x) 



(25) 



-Lp. 



For the remaining terms in Eq.(^), we can put t' — t and x' = x, and take the 
matrix trace inside the sum, 



2 cosh(f2„L) 



H{-n- 1/2) 

iJ(n+l/2) 



-iL {tanh(il„i) — S2(77i(n + 1/2)} . 

n— — oo 



We finally get 

.dG, 



(26) 



P—{P,L,fi)^-LP~iL V {tanh [{ujn - ifi) L] - sig72{n + 1/2)}. {27) 



n— — oo 



Notice that, while the first term in the r.h.s. of (|27|), coming from the singular 
behaviour of the Green function, is linear in /i, the second one, a 7r-periodic 
function of /iL, contains the finite size effects (vanishing for L — > oo). It is easy 
to see from this expression that, in the /? — s- cx) limit, the mean value of the 
particle number is 



hm N ■ 



(28) 



where [x\ means the integer part of x. 
Now, ( ^) can be integrated to obtain: 

-n G{p,L,^^)-G{P,L,{))] 



n— — OO 



f cosh UiOn ~ L] -r . r . /^x 

I coshKL] ^ +-A^^---Wr.+ l/2) 



(29) 



Taking into account that = ^'^(n-i) '^s have, for the piece of the free 
energy which depends on /i, the (finite) result 

-/3{G(/3,L,Ai)-G(AL,0)} = 



1 I 



(l + 2cos(2//L) 



(1 + 2e-2'^"-f' + e-4""-^) 
3(A*L,e-2-V/3)' 



(30) 



03(O,e-2-i/'3) 
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where 93{u,q) is the Jacobi theta function pO| . 



3 The free energy for fi = 



In order to have the complete Gibbs function, we wiU evaluate in this section 
the free energy of our system of two-dimensional fermions as a function of tem- 
perature, for ^ — 0. We take 



In Det {D{P, L)),^ = Tr { ^7°ftA;(i, t' , x') 



d_ 

dp 

where the Green function of the operator is given by (^8|) with /i = 



(31) 



fc(t, X] t' , x') — iL ^ 



(1+7^) g-la^'^H'^-^') 



2 cosh(a'„L) 



(32) 



As in Section 1, the term in the r.h.s. of Eg. (pT|) which does not contain 7^ gives 
a vanishing contribution. For the remaining, we get two contributions to the 
trace: The first, coming from the absolutely convergent series (for < x,x' < 1) 



I LXJ 



n— — 00 



H{-n- 1/2); 
2cosh(w„L) V 0; H{n + 1/2) 



L °° 

-— {tanh(tj„L) - sign{n + 1/2)} 

P n— — 00 



Op n=0 



(33) 
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The second one turns out to be divergent, and therefore a regularization proce- 
dure is required. We will use a "point-splitting" prescription to get 



-iuJnl^it — t') p — ^^OJnL{x — Xf) 



H{x' - x)H{--n - 1/2) - H{x - x')H{n + 1/2); 
0; H{x' - x)H{n + 1/2) - H{x - x')H{-n - 1/2) 

= 2^ /o dxf^ dt Im dt (h{x - x') £ e-(2"+i)f [^(--')+i/3(t-t')] 
+H{x' - x) E e"^^"+^^^ [^^'^'"'^'+'''^*"*'^] 1 

71=0 J 



|(t',a:') = (t,a:) 



sinh ( ^ [L I x - x' I - t')\ 

V VP 



> l(t',x')=(t,x) 



ttL cosh(7re//3) d ( Lf3 ttL ^ /] 

for e = L I x' — a; |<< /3, which shows a singular piece proportional to L. 
So, we have for the free energy at = 



F{(3,L) = -^ E ln{l+ e-''^"^} 
2ttL 



(34) 



= -^ln<! 



^3(0, e P ) 



1 oo 
P n=l 





ttL 


C 




+ — 
7re^ 








P 


1 7re^ 


ttL 


C 
(35) 


3/32 



with C an integration constant. 

The Casimir energy is obtained through the limit of this expression for van- 
ishing temperature, 



lim^^^F(/3,L) = EcasXL) = -i /g" da;ln{l + e'^-^} + — 



(36) 



24L 



9 



where one can see that the smgular 0(e ^) part can be ehminated through 
a renormaUzation of the vacuum energy density. The finite part of Ecas.{L), 
, coincides with the result obtained in ref . ( ) , and gives rise to an at- 

tractive force between the edges of the segment where the fermions are confined. 
Then, we get the finite result 

F{P,L)-Ecas.{L) = 



--ln<|03(O,e P ) 



P n=l 



( AttLti 
1 -e 



V 



P 24L 



(37) 

The undetermined constant C can be evaluated by imposing the vanishing of 
entropy at zero temperature. Since 



S^(3'-^F{(3,L), 
for /3 — » oo the Poisson formula for the previous series allows us to write 



(38) 



d (C 



d(3 \(3 



+ 0{p-')\ = -C + 0{(i-'), 



(39) 



implying that C — Q. 

We finally get for the free energy as a function of P,L,iJ, : 



27r f3 1 ' 



G{P,L,^i)~ EcasXL) 

ItxL 



1 oo / 4:7rLr, 

- 5] In 1 - 



ttL 



TT 

24L' 
(40) 



4 The chemical potential through a boundary 
value problem 

In this section, we will show how the result in Section 1 can alternatively be 
obtained with the methods developed in ref.|^, by introducing the chemical 
potential through suitable boundary conditions. 

As before, we will consider the differential operator 

D{P,L) = ^^°dt + ji'd,, (41) 
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with the dimensionless variables t, x taking values in the segment [0, 1]. We will 
take -D(/3, L) to act on differentiable functions satisfying the boundary conditions 

{A{r])x) {x) = x{t = l,x) + e^^xit = 0,x) = 0, 

(42) 

(Bx) [t, x) = [1 + 7i]x(i, x) = 0, for x = 0, 1, 

and call this operator {D{P, L))^^^^ g . 

The method requires a basis in the kernel of the differential operator, 
ker D{P,L). The calculation can be greatly simplified by choosing a complete 
system of functions satisfying the boundary conditions at the spatial edges. 
Such basis can be constructed from the eigenfunctions of the Hermitian opera- 
to"^ H75£)5in [0,1]: 

~*'T5^Xn(a;) = KXnix), 

(43) 

{BXn) (0) = = (BXn) (1), 

which are given by 

Xn{x) = e^'^("+V2)^75 j ^ , with A„ = 7r(n + 1/2), n€Z. (44) 

Now, a basis in ker -D(/3, L) is 

Mt,x) = e("+i/2)i[/3t+zL.75] J ^ , with neZ. (45) 

The following step is to get the projected boundary values of tpn{t,x) through 
the boundary operators ^(77) and B: 

{Ai7j)%lM) (x) \ / hn{x;v) 
Hn{t,x;r])= { {B^n){t,0) = 
(B^„)(i,l) / V 



hn{x; r}) = ipn{t = 1, x) + e"''ii)n{t = Q,x) = 

_ |^g(n+l/2)7r/3/L _|_ gj?7j gi7r(n+l/2)a:75 ( ^ 



(46) 



Forman's operator, ^(ji'^rj), is defined as 

H„{t,x;r]')=^r]',7]) Hn{t,x;r^), (47) 

for any basis in ker D[(3,L). Since the operator B does not depends on the 
parameter 77, $(77', 77) has the form 
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and our election ol basis allows us to determine ^{rj' , rj) from 

-g(n+l/2)7r/3/L _|_ ^iTjl - 



3(n+l/2)7r/3/L ^ ^iTj 



(49) 



hn{x]-q). 



Note that the hn{x;ri) ^ Xn{x), for n G Z, constitute a complete system in 
L^(0, 1), and ^{-q' ^rf) is diagonal in this basis: 



„ir,/-(n+l/2)7r/3/L 



(50) 



In what follows, we will need the quotient of $'s for two sets of parameters, 
(/3,L) and (/3o,io) respectively, 

1 _|_ g»»?'-(n+l/2)7r/3/i \ / _|_ gij;-(n+l/2)7r/3o/Lo 



($(?7 ,77) $Q (77 ,?7))„^„i - (^n.m 1 _|_ gi,,_(n+l/2)77/3/L / U + ei';'-(n+l/2)vr/3o/Lo 

(51) 

It is easy to see that [$(?]', 77) $g ^(v'tV) ~ l] is a trace class operator, so the 
determinant deti {^{r]' ,ri) $(^^(77', 77)) exists: 



deti ($(77', 77) %\il',ri))^X{ 



1 -|_ gi)7'-(n+l/2)7r/3/L 



1 + gir,-(Ti+l/2)7r/3/L 



1 



„ir;-(n+l/2)7i-/3o/Lo 



1 4- gir,/-(n+l/2)7r/3o/io 



, 7r/3 7r/3o 
93(^,e-2r) 03(§,e^2Lo) 



93(ie-2L) 03(^,e-2Lo) 



(52) 



The result established in ref.[g| is 

deti{(^(/3,L))^(^,),, (i^(/3o,io))I(,,),s (^(/^o, io))^(,),B iD{(3,L))-,\^^,,] 

deti($(r/,77) $-1(7,', r;)). 
Now, we choose the parameter rj' = ifi' P, rj = (3 and noticing that 

= (e-*''*^(/3,i)e^''*)^(o),B 



(53) 



(54) 
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we can write 



= In 



For /i' = this reduces to 

p [G{P,L,ii)- G{p,L,0)]-po [G{po,Lo,txp/Po)-G{po,Lo,0)] 



= ln. 



7r/3 . . Z§o 



(55) 



(56) 



93(^,e-2L) ^3(0,e"2io) 



and wc arc still free to choose the parameters Po and Lq conveniently. 
Notice that, for P/Po « 1, we can write 

Po [G{po,Lo,fip/Po)-GiPo,Lo,0)] = Po^p^{Po,Lo,no) 

= -11 P N{Pq,Lq,IJ,q), 

for < /Uo < /X /3//3oi and N{Po, Lq, hq) the mean value of the particle number. 
So 

Po [G{Po,Lo,fi p/Po)-G{po,Lo,0)] ^ 0. (58) 

Pq—>-oo 

On the other hand 

^^3(!iL^,e~2L^) ^ ^ 

lini 2 — ^ — ^ ^ 5g 

6*3(0,6 2Lo) 



for any value of Then, we have 

p [G{P,L,fx)- G(/3,i,0)]=ln. 



7r/3 

g3(0,e-2L) 



2vL 

-^?PL , , J ^3(0, e /5 ) 



(60) 



27r 



+ In ■ 



271-L f ' 

93(A*L,e /3 ) 
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where we have used the inversion formula for the Jacobi 03-function, 



e3{^,e-2L) = e—93ifiL,e P ). (61) 
This resuh coincides with the one obtained in Eq.(pO|). 

5 Conclusions 

In this paper, we studied the possibihty of simulating the presence of a chemical 
potential by imposing "twisted" boundary conditions in the Euclidean time 
direction. This mechanism allowed us to relate the free energy of a system of 
massless confined fermions to a functional determinant of the Dirac operator. 

To study such determinant, we made use of two approaches. The first one 
is based on the knowledge of the Green function satisfying adequate boundary 
conditions. The other relies on functional techniques Q which allow to perform 
the calculation starting from boundary values of functions in the kernel of the 
Dirac operator. 

In Section 2, the difference between Gibbs free energies with and without 
chemical potential was written in terms of the /i-intcgral of a trace involving 
the Green function. A careful treatment of this trace (which is not absolutely 
convergent), through point splitting regularization, put in evidence the presence 
of a /Lt^-dependent term. 

For completeness, in section 3, the Helmholtz free energy was computed by 
a similar method, employing this time the /^-derivative. The constant of inte- 
gration was fixed from the vanishing of entropy at zero temperature. Moreover, 
the f3 oo limit allowed us to identify the Casimir energy of the system. 

Although the results of sections 2 and 3, i.e. the Gibbs free energy of this 
simple model, could also have been obtained directly from the eigenvalues of the 
Dirac Hamiltonian (which, in this case, can be exactly evaluated), the method 
developed in those sections can be applied to more complicated and realistic 
situations, e.g. the 3 -I- 1-dimensional bag model. 

In section 4, we made an alternative calculation following the techniques 
introduced in Starting from an adequately selected basis in the kernel of 
the Dirac operator, the Fredholm determinant of quotients of Forman's oper- 
ator could be easily constructed. Through an analytic extension in the pa- 
rameter defining the twisted temporal boundary condition, we reobtained the 
/x-dependent piece of the Gibbs free energy of the system. 

The calculational scheme explored in the present paper should allow for the 
introduction of a nonvanishing chemical potential in the 3 -f 1-dimensional chiral 
bag model, thus complementing the results in |^, 0] . 
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